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Abstract 

We demonstrate the possibihty of examining cosmological signatures in the DBI 
inflation setup using the BGMPZ solution, a one-parameter family of geometries for 
the warped throat which interpolate between the Maldacena-Nuhez and Klebanov- 
Strassler solutions. The warp factor is determined numerically and subsequently 
used to calculate cosmological observables including the scalar and tensor spectral 
indices, for a sample point in the parameter space. As one moves away from the 
KS solution for the throat the warp factor is qualitatively different, which leads 
to a significant change for the observables, but also generically increases the non- 
Gaussianity of the models. We argue that the different models can potentially be 
differentiated by current and future experiments. 
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1 Introduction 

There has been much interest in recent years in cosmological applications of string theory. The 
availability of precision data relating to the cosmic microwave background (CMB) opens up the 
possibility of constraining the phenomenology of string-inspired models, in particular regarding 
the dynamics of the inflationary epoch. Essentially, signatures of the extreme conditions governing 
the behaviour of the early universe are amplifled by inflationary expansion and can be accessible 
by modern day observations. The typical energy scales involved far exceed those obtainable in 
particle collider experiments. In this paper we consider a particular inflation scenario, the Dirac- 
Born-lnfeld (DBI) model of |l|,[3|. In the usual formulation of this setup, a D3 brane rolls down a 
warped throat towards a D3 brane where the interbrane separation is identifled with the inflaton 
fleld. The brane moves relativistically, but its speed is curtailed by the warped geometry so that 
the potential energy dominates and inflation can occur. The observable consequences of this model 
depend of course upon the choice of solution for the warped throat. The original DBI proposal used 
a AdS-like geometry, with an artificially imposed cutoff where the throat joins the bulk geometry, 
which in principle is unknown although irrelevant for cosmological implications. Such a geometry is 
unstable, and more properly one should consider a solution of the supergravity fleld equations with 
non-trivial D-brane fluxes to stabilise the compact geometry by dynamically generating a cut-off. 
Not many such solutions are known analytically, but two examples are the Klebanov-Strassler 
(KS) and Maldacena-Nuhez Q (MN) results. 

Several investigations of DBI inflation from the KS geometry abound in the literature aaaaa, 
including systematic comparisons to current cosmological data. In this paper we show that it is 
also possible to examine the inflationary consequences of the BGMPZ geometries introduced in [lOl |. 
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These are a one-parameter family of geometries describing a deformation of the KS solutioqj, and 
interpolating smoothly between the KS and MN backgrounds. Although the metric is not known 
fully analytically, one can solve for it numerically and thus obtain results for numerous observables. 
The hope is that these can be compared with experimental data for quantities such as spectral 
indices, which may then allow one to extract geometrical information about the throat-geometry. 
To this aim, we show examples of how these quantities vary as a function of the BGMPZ param- 
eter ^ which characterises the family of solutions, at typical points in parameter space. We note 
that inflationary consequences of these geometries were also considered in [ill] , for slow-roll brane 
inflation rather than the DBI setup considered in this paper. 

The structure of the paper is as follows. In section 12.11 we summarise the salient details of DBI 
inflation relevant for our purposes, followed by a brief explanation of BGMPZ backgrounds. In 
section [3] we explain the details of the numerical procedure adopted in solving the BGMPZ equa- 
tions. Example results are presented in section HI and in section [5] we discuss our results before 
concluding. 



2 Theoretical Background 
2.1 DBI inflation 

The DBI inflation scenario takes place in a type IIB string theory whose ten-dimensional target 
space is compactified on a Calabi-Yau threefold with at least one warped throat, containing a D3 
brane at its tip. This then attracts a D3 brane from higher up in the throat, which rolls down the 
warped geometry. 

An approximate ansatz for the throat geometry (which can be thought of as a warped AdS space) 
is AdS^ X T^'^, with T^'^ a five-dimensional Sasaki-Einstein manifold. One may choose the metric 
to have the forrr0: 

ds^ = r^l'^{r)dx^ + f^''^{r){dr'^ + dsl), (1) 

where dx^ describes our four-dimensional space-time; r is a radial coordinate (i.e. the distance 
from the tip of the throat); /(r) the warp factor and describes the Sasaki-Einstein manifold. 
The moving D3 brane is described by its Dirac-Born-Infeld (DBI) action. To write this down, we 
use the scaled radial coordinate (p = T^/^r and warp factor /(r) = T~^f{r), where T is the string 
tension of the D3-brane. Then the action reads: 

s = -ldx' (r 1 ((1 - /02)i/2 _ 1^ _ ^(^)^ . (2) 

One now identifies (p with the inflaton field. This moves ultra-relativistically in general, but is 
subject to a speed limit (j) < f{<j))''^ arising from reality of the action. The corresponding Lorentz 
factor for the D3 brane is: 

7 = ^/l + 4M^ir^), (3) 

^The BGMPZ solutions originally arose in a different context, that of the AdS / CFT correspondence, as the 
baryonic branch of the KS solution. 

''Here we work in the Einstein frame, which is best suited to the calculation of cosmological observables. We will 
later switch to the string frame when discussing the BGMPZ solutions. 
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where Mp is the Planck mass and the Hubble parameter, with primes denoting differentiation 
with respect to (j). The fact that the speed of the brane is limited means that the dominant 
contribution to its energy comes from the potential V{(j)), and it is this that allows one to identify 
the field (j) with the inflaton. One assumes the potential to be dominated by the quadratic mass 
terrr0, and thus (following [s^) the Hubble parameter may be written: 

H{(t)) = m0(l - B^^)/Mp, (4) 

where B represents the first in a series of small correction terms arising from the kinetic energ}{§. 
It is customary to introduce the following inflationary parameters: 

_ 2M2 fH'm' _ 2Ml fH"{ct>)\ _ 2Ml (H'{ct>) i{4>)\ 



The scalar spectral index is given to linear order in these parameters by 

- 1 ~ -4eD + 27?D - 2kd- (6) 

Note that those parameters involving H"[(j)) vanish for the Hubble parameter given above with 
B = but not in the more general case where higher order kinetic and potential energy corrections 
are included. The number of e-folds before the end of inflation is given by: 

where the integral proceeds from initial to final values of the inflaton field cj). Another quantity of 
interest is the ratio of the amplitudes of tensor and scalar fluctuations, which to linear order in the 
above parameters is given by 0]: 

r:.!^. (8) 

7 

2.2 BGMPZ solutions 



The BGMPZ [10|] geometries arise as solutions in the supergravity limit of type HB string the- 
ory, in the presence of fluxes. The explicit form of the solutions for this particular one-parameter 
family has been given in [lo| in the form of solutions to a series of coupled first order ordinary 
differential equations describing the functions entering the metric tensor. The KS and MN solu- 
tions were shown to emerge as particular cases of the BGMPZ solution, which is characterised by 
a parameter ^ which smoothly in terp olates between these two extremal cases. A previously found 
deformation of the KS geometry [l5| was also shown to emerge as a particular solution. 

More specifically, the BGMPZ metric in the string frame is: 

ds^ = e^^W^fx^ dx^^ + e-^P^'^+-^'Ue + dsl (9) 



^Note that in recent explicit constructions of D3-brane inflation with embedded D7-branes such a quadratic term 
does not appear in the effective potential [l^, EI' 



^This form of the potential always fulflUs the consistency relation of [l^ . 
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with the first two factors belonging to AdS^ and ds^ describing the Sasaki-Einstein space (see ([T]) 
for the analog of this equation in the Einstein frame). It is convenient to introduce the quantity: 

v{t) = e6p(*)+2^(*) (10) 

and this together with the functions A{t), x{t) and the dilaton <I>(t) are all that is needed to examine 
the inflationary consequences of this geometry. These functions obey a set of non-linear coupled 
first-order differential equations with additional functions h2{t), a(t) arising from the fluxes and 
metric in the compact internal space. A full analytic solution of this system of equations does not 
seem to be possible, but results including appropriate integration constants are obtainable as a 
power series at asymptotically small or large t. In particular, a{t) has the form [lo| : 

a{t)^-l + Ct^ + Oit^), t^O; (11) 
^ -2e-* + at/y(-l + t)e~f +0(6"^), t ^ oo. (12) 

This defines the parameter ^ characterising the family of BGMPZ geometries. It lies in the range: 

i<«<i (13) 

with ^ = 1/6, 1/2 corresponding to the MN and KS solutions respectively. The UV {t — > oo) 
behaviour of the dilaton changes sharply as one approaches the MN solution. For ^ ^ 1/6, it 
approaches a constant, instead diverging for the MN solution itself. The parameter ajjv in equation 
(jl2p is fixed by requiring a regular solution for the metric. Similar power series are given in [Ifl] for 
the functions A{t), ^{t), v{t), h2{t) and x{t). The metric in the Einstein frame, required to make 
contact with cosmology, is given by: 

dsl = e*W/2^s2, (14) 

where $(t), as stated above, is the dilaton field and not to be confused with the inflaton field (j) 
introduced in previous sections. Combining equations ([H [9l one finds: 

y(^) ^ g-4A(t(r))-${t(r))/2^ ^^^-^ 

and: 

[■t \ ^2A{t')+x{t')+'l>{t')/2\ 2 

(16) 



r-t 

r{t)= I dt' 







v{t') 



The integration limits for r{t) are chosen to fulfil the boundary condition r(0) = 0. From equations 
([H [T5l [T6|) one sees that cosmological implications of the BGMPZ solution occur in two ways. 
Firstly, through the warp factor (including the dilaton from the frame transformation). Secondly, 
through the radial coordinate transformation r{t). The procedure for obtaining cosmological ob- 
servables is as follows. One first solves the BGMPZ equations in the string frame for the quantities 
introduced above. As discussed, this is not possible analytically and so a numerical solution must 
be performed. One then transforms the metric to the Einstein frame in order to calculate physical 
quantities. The details of the numerical solution are discussed in the following section. 
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3 Numerical Solution of the BGMPZ Equations 



3.1 String Frame Metric 

The BGMPZ equations are a set of coupled non-linear first-order differential equations in the string 
frame radial coordinate t, which are subject to a number of potential numerical instabilities around 
t = 0. Here we outline our procedure for numerical solution, before discussing the transformation 
of the resulting metric to the Einstein frame. 



We solve the BGMPZ equations [10| using a third order Adams-Bashforth algorithm. This is a 
multi-step method requiring 3 previous function evaluations at each integration step. Initial values 
are provided with a fourth order Runge-Kutta algorithm. A problem arises in that the analytic 
expressions for the derivatives of some of the metric functions are numerically unstable at t = 0. 
This is easily solved by using the power series results for a{t), v{t) etc. in this regime, and matching 
at some suitable value t = to chosen as small as possible so that the power series and full expressions 
for the derivatives are approximately equal. Furthermore, the power series for A{t), h2{t) and x{t) 
depend upon a parameter A(^) which depends on the values of the dilaton at i = 0, oo (<I>o and 
respectivelyjil: 

3(1-20 ■ 

This appears to diverge at ^ = 1/2. However, one then has the KS solution in which the dilaton 
is constant so that ^u y = ^q. Thus A(l/2) is in principle finite and is known to have a value 
A(l/2) = 0.93266 [lO|. In practise ^jjv niust be found by numerical solution of the BGMPZ 
equations, and thus equation (fT7|) becomes unstable. We solve this by interpolating the numerically 
found A over a small region close to ^ = 1/2 so that the asymptotic KS value is indeed reached. 
Numerical results for A(^) are shown in figure [1] with no such interpolation, and one sees that the 
asymptotic value is reached before the instability occurs, so that results near the KS solution will 
not be sensitive to the interpolation. For large values of t the numerical solutions for A(t) and x{t) 
also become unstable, and must be matched to the UV power series given in [l^]. The series for 
exp(2^(t)) has an overall ^ dependent normalisation which must be found numerically from the 
region over which the numerical and power series solutions overlap. In addition, some interpolation 
is necessary between the numerical and power series solutions. If An is the numerical solution and 
Ap the power series, a suitable choice is: 

e^^ ~ (1 - 0.5e"(*-*^))e2^" + 0.5e°(*-*l)e2^^ t < h; 

~ 0.5e-"(*-*i)e2^" + (1 - 0.5e-°(*-*i))e2^p, t > h, (18) 

with a selected to ensure a smooth interpolation. We have checked that final results for those 
quantities depending on the derivative of the warp factor are not sensitive to this interpolation. 
A further check on our numerical solutions can be performed by comparing to the known KS and 
MN solutions at ^ = 1/2 and ^ = 1/6 respectively. We have explicitly verified that the shape and 
normalisation of the analytical results for the above functions in these limits are reproduced. 



^One is free to choose the initial value of the dilaton, as in the supergravity approximation this amounts merely 
to an overall scaling of the metric. We set "I>o = for our analysis. 
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Figure 1: Behaviour of the parameter A(i^) defined in equation (|17p . together with its asymptotic 
value A(f/2) = 0.93266. One sees that this value is reached before the numerical instability sets in. 



3.2 Transformation to the Einstein Frame 

Once the BGMPZ metric has been found, equations (jl51 [TBI) provide numerical solutions for the 
Einstein frame warp factor and radial coordinate. The numerical solution for these functions con- 
sists of a series of values at discrete values of t. Values at general values of t are obtained by 
interpolating between those on the discretised axis of t points, and a linear interpolation is found 
to be sufficient. Equation (fT6]) must be integrated numerically, and we use a 96-point Gaussian 
integration. 

When calculating cosmological observables, one must find the value of the string frame coordinate 
t that corresponds to a given value oi 4> oc r(t). Thus the solution for r{t) must be inverted numer- 
ically to find t{r). We use a bisection algorithm with sufficiently many iterations. 

Each of the numerical procedures described above introduces various tuning parameters (e.g. step 
sizes, interpolation parameters, iteration numbers) which control the convergence and accuracy 
of the final numerical results. We have checked the (in)sensitivity of all final results to these 
parameters. 
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4 Results 



4.1 Warp Factors 

Several parameters must be fixed in order to produce numerical predictions from the above analysis. 
Firstly, the values of the string scale and inflaton mass. Secondly, there is an overall normalisation 
of the warp factor /(</>), which as well as the string scale and coupling has a dependence on the 
D3 brane tension and the deformation parameter of the conifold solution. Furthermore, there is a 
scaling factor TZ in the definition of the inflaton (j), which also depends on the conifold deformation, 
brane tension and, in the case of inflation with higher dimensional wrapped branes , upon 

a scale factor involving the volume of the wrapped cycle. 



Here we present results for the following choice of parameters: 

nis = O.OlMp; m = IQ-^Mp- /(O) = 10^ 7^ = 25, (19) 

where is the string scale (we work in units such that = 1), m the inflaton mass and 
<j){t) = TZr(t) . Our objective here is merely to demonstrate an example of how the numerical 
solution of the BGMPZ equations can be applied to calculate cosmological observables. We discuss 
in more detail the possible constraints on DBI models in section 14.41 

The BGMPZ warp factors for ^ = 1/2, 0.167 (corresponding to the KS solutions and a solution 
close to the MN case respectively) are shown in figure [2] as a function of (j)/(f)e, where (pe is the value 
of the inflaton field at which the warped throat joins onto the bulk geometry defined by: 

fi<Pe) = 1. (20) 

Thus, (j)e has a different value for each of the geometries. One sees this from figure[3l which shows (pe 
as a function of ^. One sees that (j)e diverges towards the MN solution (^ — > 1/6) as expected from 
the divergence of the dilaton $(t). Two things are evident from flgure[2j Firstly, both warp factors 
converge for low values of </), corresponding to low values of the string frame radial coordinate t 
where the geometries are the same. Constancy of the warp factor as — > demonstrates that the 
BGMPZ solutions do indeed have a dynamically generated cutoff of the warped throat. Secondly, 
the warp factors converge as (j)/(j)e 1- This is a consequence of both geometries describing warped 
throats that are asymptotically AdS. 



The difference between the warp factors in figure EJ which effectively represent extremal cases of 
the BGMPZ geometry, emerges in the shape at intermediate (p values. The KS throat (^ = 1/2) 
interpolates smoothly between the fiat warp factor at low values of (p, and the asymptotic AdS 
regime towards the tip. It can be effectively parameterised by a function of form [H]: 

f(4>) = — 4> = ^, (21) 

(;u2+,/.2)2' Cj)e 

for suitable parameters Af and /i. An extra qualitative regime is present in the general BGMPZ 
warp factor, namely a shoulder at intermediate (p values whose slope differs from the eventual 
asymptotic AdS behaviour. Thus, an extra parameter is needed in the denominator of equation 
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Figure 2: Warp factors from the numerical solutions of the BGMPZ equations. Results are shown 
for ^ = 1/2 (solid line) and ^ = 0.167 (dashed line). 
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Figure 3: Value of the inflaton (/)e corresponding to the radial coordinate at which the throat is 
glued to the bulk geometry, defined by f{4>e) = 1- 
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I 

Figure 4: Visualisation of the geometry of two throats at (a) ^ = 1/2 and (b) ^ = 0.167. The 
horizontal distance is in units of </>/0e) and the proper distance between two points at different 
values of the inflaton field is represented by the length of a geodesic path connecting them along 
the throat surface. 

(pTI) , and an appropriate generalisation of the function (f2T]) was also first suggested in : 

f((j)) = i ^. (22) 

We have checked explicitly that this function is able to reproduce well the shape of the BGMPZ 
warp factor. Thus, the results of this paper confirm that phenomenological warp factors of form 
()22|) do indeed correspond to geometries exactly realisable in string theory - namely the BGMPZ 
solutions. 

To give an impression of how one might visualise the two extremal cases at ^ = 1/2 (KS) and 
^ = 0.167 (close to the MN solution), Figure [4] shows two-dimensional slices of the geometry, 
parametrised by the radial parameter (p/4>e- The vertical direction is used to show the warping, in 
such a way that the proper distance between two points at different radial values (given by equation 
([1])) is represented by the length of a geodesic path connecting them on the throat surface. The 
heights of the throats are scaled by (pg for each case. The KS throat is then longer on such a plot 
than the MN-type solution, consistent with figure [21 due to the higher warp factor. 

Although viewing the warp factor as a function of (f)/(f)e (as in figures [2] and H]) allows one to 
more easily visualise the geometries, the physical scale on which to consider the warp factor is the 
number of e-folds. Thus, in figure [5] we show the warp factors for ^ = 1/2, 0.167 as a function of 
Ng. Contrary to the naive expectation of figure [21 the warp factor at a given value of A'^g increases 
as one moves away from the KS solution. This is due to an increase in 0e (and, thus, the value of 
(f) at which a given N^. occurs) as ^ decreases. 
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Figure 5: Warp factors from the numerical solutions of the BGMPZ equations considered as a 
function of the number of e-folds before inflation ends, N^.. Results are shown for ^ = 1/2 (solid 
line) and ^ = 0.167 (dashed line). 
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Figure 6: Scalar spectral indices from the numerical solutions of the BGMPZ equations, followed 
by transformation to the Einstein frame. Results are shown for ^ = 1/2 (solid line) and ^ = 0.167 
(dashed line). Also shown are la and 2a error bands from the WMAP3 data. 



4.2 Scalar Spectral Index 

In figure [6] we show the spectral index as a function of the number of e-folds before inflation ends, 
defined by equation ([7]). Noticeably, the qualitative behaviour of the spectral index at high e-fold 
numbers differs between the two cases. The KS solution shows a bump which is absent as ^ tends 
towards the MN background. Furthermore, this bump occurs (for our choice of parameters) within 
the last 50 e-folds of inflation. Superimposed on the plot are la and 2a error bands for the spectral 
index Ug at 0.002/Mpc as measured from current data including WMAP3 [l8| : 

= 0.947 ±0.015. (23) 

Note that this value is derived from a fit to combined data sets performed under the assumption 
that the ratio of tensor to scalar fluctuations, r, is zero. More properly one should should consider 
confidence contours in the (r, n^) plane. Then, for non-zero r, the central value of Ug increases 
but the uncertainty band has roughly the same width for reasonable variation of (see figure 14 
in [3]). In the DBI scenario, r is indeed non-zero and varies with ^. Nevertheless, we plot the 
value of equation (|23|) in figure E] so as to indicate the amount of uncertainty encountered in current 
measurements of n^. The uncertainty band in any case moves relative to the curves as a result of 
varying the DBI parameter space. 

In figure [7] we show the value of jn^ — 1| at A'^e = 55, as a function of the ^ parameter characterising 
the BGMPZ geometries. For the = 0.167 case, one sees two ranges of in which the value of 
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Figure 7: Scalar spectral index 55 e-folds before the end of inflation, as a function of the BGMPZ 
parameter ^. 



\ns — 1| is less than and greater than the KS result respectively. As ^ increases from 0.167, the 
crossing point of the two curves also varies such that for low ^ it lies at A^e < 55 whereas for higher 
it lies at > 55. Thus, for higher ^ the value of In^ — 1| increases with increasing ^, whereas 
the opposite is true at low ^. For points in the parameter space other than those of equation (fT9]) . 
the spectral index may either increase or decrease as one moves away from the KS solution. 

It is also interesting to note from figure [7] that the value of |ns — 1| at Ne = 55 for the MN-like 
solution lies outside the la WMAP3 error band such that is slightly disfavoured relative to the KS 
solution. For other points in the parameter space, this difference may be enhanced and thus the 
ability of current data to constrain the ^ parameter along with the other parameters is not in doubtH. 

Also of interest is the running of the spectral index. However, for the parameter space point chosen 
for our illustrative results, the running is sufficiently small for each geometry as to be unobservable. 
A plot of dus/dln k for the KS solution can be seen, albeit for a different region of parameter space, 
in H- 

*We found examples of other points in parameter space where the KS result for \ns — Ij lies outside the 2a error 
band from WMAP3 whereas the MN result does not, although these choices of parameters were also ruled out by 
non-Gaussianity constraints, discussed in section 14.41 
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Figure 8: Ratio of tensor to scalar modes from the numerical solutions of the BGMPZ equations, 
followed by transformation to the Einstein frame. Results are shown for ^ = 1/2 (solid line) and 
^ = 0.167 (dashed line). 



4.3 Ratio of tensor to scalar fluctuations 

One may also solve for the ratio of tensor to scalar modes r, given in equation Its behaviour 
as a function of the number of e-folds before the end of inflation is shown in figure [8] for = 1/2, 
0.167 where analogously to the case of the spectral index one sees a markedly different qualitative 
behaviour. The value of r at A'^g = 55 is shown in figure [9] as a function of Similarly to the 
spectral index, one sees that r may increase or decrease as ^ increases. 



Here we only wish to illustrate the sort of variation that can arise from the variation of the BGMPZ 
parameter, and the values of r shown in figures [8] and [9] cannot represent predictions without first 
ensuring that all constraints on DEI inflation are obeyed (see the following section for a discussion, 
and [l^ for a discussion of gravitational waves in DBI inflation). A full imposition of all theoretical 
and experimental constraints is beyond the scope of this paper. Nevertheless, in figure [10] we 
show the value of r versus Ug at 55 e-folds before inflation, for various values of ^ and using the 
parameters of equation (I19p . Comparing with the corresponding figure in [18l], one again reaches 
the conclusion that experimental data can constrain the BGMPZ parameter taken alongside the 
rest of the DBI parameter space. Theoretical constraints would act to further reduce the available 
parameter space. 
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Figure 9: Ratio of tensor to scalar modes 55 e-folds before the end of inflation, as a function of the 
BGMPZ parameter ^. 
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Figure 10: Ratio of tensor to scalar modes vs. the scalar spectral index at 55 e-folds before inflation 
for (a) ^ = 0.167; (b) ^ = 1/2. Points on the curve are spaced equally in ^. 
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4.4 Constraints on Parameters in DBI Inflation 



In this paper we have presented results for a particular point in the parameter space of the DBI 
inflation model, namely that of equation (jl9p . This is merely to show that it is possible to calculate 
cosmological observables for this class of throat geometries, and also the possibility that geometries 
with different ^ values could in principle be differentiated by current and future data. A much more 
thorough analysis would involve a scan over the parameter space including the various theoretical 
constraints that are known to affect DBI inflation models, as it has been done for a generic form 
of the warping in [9(]. 

There are a number of theoretical constraints which act on the parameter space of a given DBI infla- 
tion model. Firstly, there is the requirement that the supergravity approximation used in deriving 
the geometry is valid, which includes the condition that one may safely ignore the back-reaction 
of the colliding D-branes on the throat. There is also the requirement that the throat fit inside 
the bulk geometry. More practically, one must ensure that sufficiently many e-folds of inflation are 
generated. 



A further concern with DBI inflation in particular is that it generates a high degree of non- 
Gaussianity in the power spectrum of cosmological perturbations, as a direct consequence of the 
large Lorentz factor associated with ultra-relativistic motion of the moving D-brane 0, [sS] • It has 
been argued, for example, in 

@, E HH that present non-Gaussianity measurements can be used 
to heavily constrain the parameter space of the D3-D3 model by requiring ''){N(. = 55) < 31. This 
is not a problem for the parameter choice of equation (jl9p , as seen in figure [11] where we show the 
behaviour of j{N(.) for = 0.5 and ^ = 0.167. One sees that both the KS and MN-like results lie 
well within the upper bound from non-Gaussianity, but that 7(0) rises sharply as the number of 
e-folds decreases such that this constraint is a significant one. We also note that 7((/>) at a given 
A'^e increases as ^ decreases, in other words that the non-Gaussianity increases as one moves away 
from the KS throat. 



The non-Gaussianity constraint can be relaxed by considering generalisations of the original DBI 
scenario. One such extension is to replace the D3 brane with higher dimensional D-branes that 
wrap cycles in the throat , [13] , which might relax the previous constraints on the DBI inflation 
parameter space. In principle, extended D-brane scenarios could be embedded in the family of 
throat geometries discussed in this paper, although the question of how large an effect the back- 
reaction of higher dimensional branes has on the inflaton potential or in constraining the DBI 
parameter space remains uncleaijfl. Regarding non-Gaussianity, it is interesting that moving away 
from the KS throat results in a 7 factor which is higher for a given number of e-folds. This suggests 
that analyses of DBI inflation based on the Klebanov-Strassler throat are underestimating possible 
non-Gaussianities. 



Another influence on the parameters comes from the possibility that one could deal with multi- 
field infiation. In this work we have only taken one direction of the D brane as the scalar inflaton 
field, but this could be extended to include more directions. As discussed recently in 0, HI], it can 
transpire that multi-field models are well-approximated by single-field limits, which can also have 

'a preliminary study of back-reaction effects in higher dimensional brane inflation scenarios is presented in [l^ . 
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Figure 11: Behaviour of the Lorentz factor 'y{(j)) as varies, for ^ = 1/2 (sohd) and = 0.167 
(dashed). Also shown is the upper bound jiNe ~ 55) < 31 from non-Gaussianity constraints. 



an effect in relaxing previously given constraints on DBI models. 

One has to be careful however assuming that it is possible to vary the parameter ^ completely 
independently. A priori the potential will change as we deform the geometry and therefore the 
ranges of all other parameters will be influenced. In our simplified approach, we did not take these 
effects into account. 

In light of the uncertainties above, the example parameters (|19p used in this paper may or may 
not be ruled out by present measurements. However, we seek to demonstrate only the fact that 
cosmological observables do depend quantitatively on the deformation parameter ^ of the BGMPZ 
geometries, and it is still very likely to be the case that generalised DBI setups can be approximated 
to a reasonable degree by a throat of this kind. 



5 Discussion 

In this paper we have explicitly demonstrated the possibility of calculating cosmological observables 
in the DBI inflation setup using a one parameter family of type IIB supergravity solutions that 
describe the geometry of a warped throat, the BGMPZ solutions of Butti et al. ^IQ] that interpolate 
smoothly between the Klebanov-Strassler and the Maldacena-Nuiiez solution. We have provided 
examples of cosmological parameters, namely spectral indices, that can be calculated from the un- 
derlying geometry. 
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The solution for the metric of the geometries in question is not possible analytically. Therefore 
numerical methods have been used and shown to provide an adequate representation of the metric 
in terms of numerical precision. Instabilities in the derivatives entering the equations have been 
dealt with by matching to known power series expansions of the metric functions at asymptotically 
small and large values of the radial coordinate. 

We presented warp factors for two almost extremal ends of the family of solutions parametrised 
by ^: the Klebanov-Strassler throat = 1/2), and a geometry close to the Maldacena-Nuiiez 
solution = 0.167). The qualitative behaviour of the warp factors is seen to be different. Both 
solutions show a flat warp factor at low values of the rescaled radial coordinate (/>, corresponding 
to a dynamically generated cutoff, and an asymptotically AdS behaviour at large values as ex- 
pected. However, the warp factor as it moves away in ^ from the KS solution develops a shoulder 
at intermediate values, whose slope is in general different from the eventual asymptotic behaviour. 
This can be effectively parameterised by a function of form (I22p . thus explicitly demonstrating 
that phenomenological fits of DBI inflation to cosmological data that assume a warp factor of this 
form indeed correspond to exactly realisable warped throat geometries in a known string theory 
compactification. 

We presented examples of scalar spectral indices and the ratio of tensor to scalar modes calculated 
using the different geometries. The different qualitative behaviour observed in the warp factors 
carries through to the spectral indices, and a quantitative difference between the values of the 
indices from different solutions at 55 e-folds before the end of inflation is potentially measurable. 
Constraints on non-Gaussianity and measured values of the scalar spectral index could already be 
used to rule out some regions of ^, when it is considered alongside the rest of the DBI parameter 
space. Very generically, we have found that the amount of non-Gaussianity increases as one moves 
away from the KS solution (all other parameters being equal). 

Since the warping changes qualitatively as one varies ^, it is also conceivable that that certain 
constraints on the parameter space may be relaxed by including ^ as an additional parameter with 



respect to existing analyses such as |23|. This issue certainly deserves further study. 
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